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ABSTRACT 


* 


This  paper  discusses  the  classical  Lamb  problem  for  the  elastic  wave 
equation.  The  motivation,  for  the  authors,  is  to  be  able  to  conveniently 
construct  Green's  functions  (matrices)  for  later  use  in  formulating  and 
solving  various  inverse  problems.  For  example,  we  will  want  to  be  able  to 
solve  for  perturbations  from  constant  reference  densities  and  LamA 
parameters.  Hence,  Green's  function  for  the  homogeneous  isotropic  equation 
is  discussed.  Due  to  the  scattering  taking  place  in  inverse  problems  it  is 
usually  impossible  to  retain  the  P-SV  and  SH  decoupling;  hence,  we  do  not 
pursue  this  decoupling  in  the  formation  of  the  Green's  functions  herein. 

While  nothing  conceptually  new  is  presented  here,  the  approach  is  a  bit 
different  and,  we  believe,  is  helpful  in  isolating  some  important  issues. 
The  approach  is  algebraic  in  nature  and  makes  heavy  use  of  several  simple 
facts  from  linear  algebra;  for  example,  the  spectral  decomposition  of 
special  matrices.  This  approach  facilitates  some  helpful  decoupling, 
particularly  in  solving  for  reflection  coefficients. 


GLOSSARY 


a,  b ,  c 


A 

B(  • ) 

cp*  cs 
C 


D.  D' 


f 

F 


H 

I 

k 


k 

A 

k 


k' 

L(  • ) 


generic  vectors  in  R* . 

(varions)  3  by  3  matrices. 

boundary  operator  at  x,  =  0  surface;  equation  (17). 
pressure  and  shear  speeds;  above  equation  (12). 

3  by  3  matrix  used  in  constructing  R;  equation  (28). 
determinants  associated  with  reflection  coefficients; 
equation  (36). 

constant  body  force;  equation  (2). 

free  space  portion  of  Lamb  Green's  function;  equation 
(18). 

arbitrary  body  force;  equation  (1). 

free  space  Green's  function;  equation  (11). 

portions  of  G;  equation  (12). 

vector  associated  with  horizontal  shear  wave; 
equation  (8b). 

Green's  function  for  Lamb  problem;  equation  (37). 

3  by  3  identity  matrix;  equation  (5). 

=  (k1,k,,kI)  spatial  Fourier  transform  variables; 
equation  (5). 

-  |k |  ;  below  equation  (5). 

=  k/k  ;  below  equation  (5). 

=  (k1,k2)  ;  equation  (18). 

basic  differential  operator  in  matrix  fcrm;  equation 
(5). 


-  i  i  - 


M(  • ) 


N».  N, 
2 


operator  L(-)  -  pu>*I  ;  equation  (9). 
portions  of  M-1;  below  eqnation  (11). 
vector  associated  with  pressure  wave;  eqnation  (8a). 
reflection  coefficients,  subscripted  with  p,  h,  v; 
equation  (32). 

reflection  portion  of  the  Lamb  Green's  function; 
equation  (28). 

transpose  of  vector  or  matrix, 
displacement  vector;  equation  (1). 

vector  associated  with  vertical  shear  wave;  equation 
(8c). 

=  (li.ij.i,)  spatial  coordinate;  equation  (1). 

=  (x1(i,);  equation  (18). 

Dirac  delta  function;  equation  (2). 
gradient  operator;  equation  (2). 

Laplacian;  equation  (3). 
small  positive  number. 

Lamd  parameters;  equation  (1). 
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density;  equation  (1). 
eigenvalues;  equation  (7) 
frequency:  equation  (2). 
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INTRODUCTION 


The  purpose  here  is  to  present  a  somewhat  different  point  of  view 

toward  constructing  Green's  functions  (matrices)  for  a  variety  of  problems 
associated  with  the  elastic  wave  equation.  The  primary  interest  of  the 

authors  is  inverse  problems  for  various  acoustics  and  elastic  experiments 
(e.g.  Cohen  et  al  [4]  and  Boyse  and  Keller  [2]).  In  this  regard  a  clear 

view  of  certain  basic  principles  and  a  systematic  approach  to  constructing 

Green's  functions  is  particularly  important.  A  convenient  setting  for  such 
studies  is  the  frequency  domain,  hence  the  emphasis  below  is  in  obtaining 
the  corresponding  Green's  functions.  However,  some  mention  is  made 
regarding  techniques  for  inverting  to  the  time  domain. 

The  approach  taken  here  is  algebraic  in  nature,  making  repeated  use  of 
several  very  elementary  properties  of  square  matrices.  Consequently,  it  is 
a  bit  more  contemporary  in  style  than  that  found  in  most  sources  and,  we 
believe,  a  little  more  natural.  There  are  of  course  many  high  quality 
studies  of  the  classical  problems  discussed  here  (e.g.  Aki  and  Richards  [1] , 
Cerveny  and  Ravindra  [3]  and  Johnson  [7])  and  the  authors  lay  no  claim  to 
the  discovery  of  heretofore  unknown  phenomena,  or  even  methods.  The  hope  is 
that  the  general  approach  will  be  useful  in  attacking  new  problems  in  a  more 
straight-forward  fashion. 

Rather  than  attempt  to  discuss  things  in  great  generality,  we  choose  to 
illustrate  the  approach  on  a  well-known  problem.  In  particular,  our  goal 
here  will  be  the  construction  of  the  Green's  function  for  (one  statement  of) 
the  Lamb  problemt  a  single  traction-free  reflecting  surface  at  z  *  0. 
Along  the  way  the  "free  space"  Green's  function  will  be  obtained.  In  order 
to  exploit  the  convenient  algebraic  properties  we  tend  to  Fourier  transform 


as  much  as  possible  at  each  stage.  By  making  use  of  a  simple  algebraic 
decomposition  of  a  constant  matrix  (associated  with  P.  SV  and  SH  modes) 


IV 


<'* 

(V 

1 

R 


valuable  insight  is  gained  and  a  certain  amount  of  desirable  decoupling  is 
facilitated.  For  example,  the  nine  reflection  coefficients  are  rather 
naturally  found  by  solving  three  3  by  3  linear  systems  (instead  of  one  9  by 
9  system) . 


SOMB  PRELIMINARIES 


The  linear  equations  of  elasticity  in  an  isotropic,  inhomogeneous  media 
are  given  in  terms  of  the  displacement  vector  u  =  u(x,t)  as 

p3*  u=  V(XV-u)  +  V  X  (MV  X  u)  +  2(V-MV)u  +  g(x.t)  , 

(1) 

x  =  (xi,xj,xj)  ,  t  >  0 

In  particnlar,  we  seek  various  solutions  when  g(x,t)  =  8(x-y)8(t)f  where  f 
is  a  constant  vector  or,  later,  a  constant  matrix.  Bence,  we  seek  various 
Green's  functions  for  this  equation.  We  Fourier  transform  in  time 
(t  ->w,  o(x,t)  ->  u ( x , (i>) )  producing 

V(XV-u)  +  V  X  (MV  X  u)  +  2 (V*  pV)u  +  p<A  =  -8(x-y)f  .  (2) 

Our  attention  will  focus  on  the  homogeneous,  isotropic  case  in  which  X, 
M  and  p  are  constants.  Then  (2)  becomes 

(X+p)V(V*u)  +  pAu  +  pw*u  =  -6(x-y)f  .  (3) 


We  use  the  following  notation: 

1.  u,  x,  etc.  will  denote  column  vectors.  When  a  row  vector  is 

needed,  we  write  e.g.  uT.  where  T  denotes  transpose. 

T 

2.  The  inner  product  of  two  vectors  is  denoted  by  uAy  or  u*v.  The 

outer  (or  tensor)  product  is  denoted  by  uv1  i*  of  course  a  3 

by  3  matrix. 

3.  k,  etc.  will  denote  unit  vectors*  k  =  k/|kj>  Moreover,  k  denotes 

Ifcl- 


-  3  - 


,sr. 


4. 


k'.  x’,  etc.  will  denote  the  first  two  components  of  k.  x; 
e.g.  k'  =  (k^k,). 

5.  Fourier  transform  is  defined  by 

f(t)  ->J  dt  el6)t  f(t)  . 

—  00 

00 

g(x)  — dx  e  -  g(x)  ,  x  e  R*  or  R* 

—  CD 

(Note  the  difference  in  definition  for  time  and  spatial 
transforms.  This  facilitates  the  notion  of  "outgoing"  below.) 

Some  elementary  facts  from  linear  algebra  are  now  stated.  These  facts 
are  established  in  about  any  standard  text  (e.g.  Strang  [9],  Hoffman  and 
Kunze  [6],  or  Mostow  and  Sampson  [8]),  and  can  be  easily  verified  by  the 
reader. 

1.  Let  a  be  a  unit  vector.  The  matrix  P  =  aa^  is  a  (orthogonal) 
projection  matrix  (i.e.  P*x  =  P(Px)  =  Px;  moreover,  if  x  is 
orthogonal  to  a  then  Px  =  0).  Clearly  Px  =  aa^x  =  ca ,  c  =  a*x;  so 
P  projects  R1  orthogonally  onto  the  line  spanned  by  a. 

2.  Suppose  a,  b  and  c  are  unit  vectors  and  mutually  orthogonal.  Then 

the  identity  matrix  can  be  decomposed  into  the  three  projections; 
T  T  T 

I  =  »»  +  bb  +  cc  .  (This  is  easily  verified  by  applying  both 

sides  to  an  arbitrary  vector  x)  . 

3.  Suppose  A  is  a  3  by  3  matrix  with  orthonormal  eigenvectors  a,  b .  c 
and  eigenvalues  o1 ,  o,,  o, .  respectively.  Then  A  has  the  spectral 
decomposition 


T  T  T 

A  =  a  aa  +  a  bb  +  a  cc 

1 -  1 -  3 - 


(This  is  verified  by  applying  both  sides  to  an  arbitrary  vector 


x  =  ata  +  a,b  +  a3c) .  Finally,  if  the  oj  *  0  then  A  is  invertible 
and  A  1  has  the  convenient  expression 


,-i  1  T  1  T  T 

A  =  —  aa  +  —  bb  +  —  cc 


a  —  a  —  <y  — 
i  a  3 


(This  is  verified  by  multiplying  AA  1  and  using  fact  2.  above) 


4.  For  vectors  in  Rs ,  the  following  is  easily  verified: 


T  T  T  T 

[u.v.w]  b  =  ua  +  vb  +  wc 


5.  Suppose  Ab  =  <rb ,  then  observe  that 


(A  +  el)b  =  (<r+e)b 


i.e.  the  eigenvalues  of  A  +  el  are  those  of  A  shifted  by 


e,  and  the  eigenvectors  are  unchanged. 


Returning  to  the  elastic  equations  (3),  we  next  Fourier  transform  the 


spatial  variables  (x  ->k)  thus  reduc ing  ■'sipr  problem  to  an  algebraic  one  (at 


least  temporarily).  The  resulting  equation  can  be  written 


(L(k)-p<.)*I)  u(k ,w)  =  e_1-'y 


where 


m 


L(k)  =  (X+p)  kk  +  jjk  I 


2  »  <T  2 

=  (X+p)  k  kk1  +  pk  I 
=  A  +  pk*I  . 

Formally,  therefore,  if  no  boundary  conditions  are  imposed,  the  ’’free" 
solution  in  the  transform  variables  can  be  expressed 

u(k.u)  =  (U^-pc/l)-1  e~  i-  ’  ^f  (6) 

In  the  next  section  we  invert  this  (k  -> x )  in  the  simplist  possible  context, 
obtaining  the  "free  space"  Green's  function.  This  will  make  good  use  of  the 
above  algebraic  observations  and  will  set  the  stage  for  the  subsequent 
section  where  the  Lamb  problem  is  considered. 


-  6 
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THE  FREE  SPACE  GREEN'S  FUNCTION 


By  the  free  space  Green's  function  we  mean  the  matrix  solution  to  (3) 


1 

{ 


with  f  =  I.  Perhaps  the  easiest  way  to  invert  (k  — >  x)  the  expression  for 
u(k,(o)  in  (6)  is  to  first  perform  the  spectral  decomposition  of  the  matrix 
in  (5).  This  turns  out  to  be  surprisingly  easy  and  generally  informative. 

First,  note  that  it  follows  immediately  from  the  facts  noted  above  that 
the  matrix  A  =  L(k)-nk*I  =  (X+p)kk^  above  is  a  projection  operator  (onto  k): 
hence  k  is  an  eigenvector  with  associated  eigenvalue  (A.+p)k*.  Moreover,  the 
only  other  eigenvalue  is  zero  and  the  associated  (2-dim)  eigenspace  must  be 
orthogonal  to  k.  It  follows  that  L(k)  =  A  +  pk*I  has  the  same  eigenvectors 
and  (shifted)  eigenvalues 


cri  =  (X+p)k2  +  pk*  =  (A+2p)k1  , 


<ri  =  jik 


For  convenience  we  select  and  denote  the  eigenvectors  as  follows. 
Associated  with  ax  we  take 


* 

2  =  *  =  k/  |k| 


where  p  suggests  pressure  (or  longitudinal)  wave.  For  a,  can  choose  any 
two  vectors  orthogonal  to  j>.  We  select 


•a 

#1 

!sl 


k*  +  k* 


(k  ,-k  .0) 
l  a 


7T  (k,.-kl,o) 


v  =  —^7  (k  k  ,k  k  ,-lk'  I*)  . 

-  VV '  II  21  1  I 


‘/.•.‘i  \*  v'. 


M  | 


Note  that  p,  ^  -  v  form  as  orthonormal  set.  The  symbol  h  suggests  the 

horizontal  component  of  the  shear  (or  transverse)  wave  and  v  suggests  the 
vertical  component  of  shear.  In  particular,  h  is  horizontal  in  the  sense 
•(0,0,1)  =  0.  In  the  homogeneous  case  (p,  X  and  p  constant  and  no  body 

forces),  the  corresponding  waves  will  remain  orthogonal  until  encountering 
some  form  of  inhomogeniety  or  boundary. 

We  can  now  decompose  the  matrix  L(k)-pw*I  using  these  eigenvectors  and 
eigenvalues  as  follows 

M(k)  =  M(k ; u>)  «  L(k)  -  pw*I 

=  (o^  -  pw  )  pp1  +  (oj  -  pw  )  Ihh  +  vv  I 

Note,  once  again,  the  eigenvalues  for  M  have  been  shifted  by  pw*  relative  to 
those  of  L.  Assuming  -  pw*  *  0  and  a2  -  pw*  ^  0,  the  inverse  of  M(k)  is 

M(k)  1  =  - - — -  ppT  +  - - — —  FhhT  +  vvT]  .  (9) 

o^-pw  9}-p(i)  L 

Hence,  the  formal  "free"  solution  to  (3)  (in  (k,w)  -  space)  is 

u(k.w)  =  M(k)-*  e"1-’*  f  (10) 

where  M(k)-1  is  given  by  (9).  To  complete  the  construction  of  the  free 
space  Green's  matrix  G(x,t;y)  it  remains  to  invert  transform  the  matrix  in 
( 10) ,  i . e .  matrix 

G(k,<o;y)*e*-^M(k)  .  (ID 

This  inversion  follows  familiar  procedures  (e.g.  see  Aki  and  Richards  [1]), 


fl 

'iS 

W 


but  is  outlined  here  for  sake  of  completeness . 

In  (9)  we  replace  hh^  +  vy^  by  I  -  pp^  (algebra  fact  2.  above)  and 
using  (7)  we  have 


(k)  1  =  — L_  ppT  +  — L__  fi  _  eetJ 

n  -nw  ft  -ftw  *■  J 


O^-pU)  O^-pU) 


1  „  [  1  1  T 

=  - -  j  +  - - - -  pp 

Cj-pto  L  Oj-pu  o4-pii> 


a  a 


i  if  c  -  c 

— -  I  +  -  - ? - 2 -  ki 

:W>  P  (cVVHcVV)  ~ 


=  Nt(k)  +  N^k) 


where 


c*  =  hlhL 

p  p 


a  u 
C  =  . 


s  p 


Hence  in  this  notation. 


G(k.u>;y)  =  e  -Z  [^(k)  +  N^Ck)] 


s  G^k.wjjr)  +  Ga(k,w;j) 

The  inversion  (k  ->x)  of  Gx  is  standard  and  proceeds  or  follows: 
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’«,**»  T  --  «  -  *  '  N.  -  •  1  . 

,0  ,*•  /V*v\\V 


G^x.w.y) 


1  f  ik-x  -ik-y  1  T 

77 — 7  dk  e  -  -  e - n - —  I 

(2n)  P(csk  "  w  > 


- i -  f  dk  e  "  (-  -)  _ I -  I 

Unj’pc1  J  k*  -  «  /c* 

s  -<*>  s 


This  triple  integral  can  be  evaluated  by  going  to  spherical  coordinates 
k  — Mr, 0,p).  These  steps  are  outlined  below  in  discussing  Ga.  The  result, 
after  using  p  =  pc|.  is 


Vz'w;z>  =  4^rFFTexp  [f  1 


The  more  interesting  step  is  the  inversion  of  Gj  in  (12),  which  we 
outline. 


G  (x.w;y)  =  — —  [  dk  e*-  ^ 

*  "  '  (2n)  J 


2  2 
C  -  C 

- 5 - 5 - ki 

p(cY-  w*)(cV-  w1) 

P  S 


c*  -  c*  f 

=  JL  _ E 

8n*  p  c*c*  J 


-£ _ E  dk  e^'D  _ I -  kk1 

p  cV  J  "  (k*  -  (0*/c*)(k*-  w*/c  *>" 

s  p  -«  P  S 


a  a 

c  -  c 
1  s  p 

a  i  a  a 
8n  pc  c 


)*  f  e**” (i"y) 

[•*7  J  “  <k*  -  .*/«*)(»■  -  .Vi 


Note  that  the  term  in  brackets  is  a  matrix,  say  B  =  (bij),  with  elements 


iHM 


-V v.V  ■  V 


i  j  J  (k  -  u>  /c  )  (k  -  w  /c  ) 

-®  p  s 


Observe  that  this  integral  is  like  the  integral  for  G  in  (13).  except  for 
the  additional  poles  at  k  =  ±  w/cs.  In  evaluating  the  integrals  in  (13)  and 
(IS)  the  usnal  procedure  is  as  follows.  Con\ art  to  spherical  coordinates, 
k  — Mr, 0,0),  with  r  =  III  and  x-jr  serving  as  the  "north  pole."  The  angular 
integrals  are  done  routinely;  then  the  r  integral  is  done  by  the  residue 
theorem.  As  in  Aki  and  Richards  [1],  one  assumes  Im(o>)  2  0  which  dictates 
that  one  take  the  residue  integration  path  with  Im(r)  2  0.  This  results  in 
the  two  poles  at  k  =  +  w/cp  and  k  =  +  w/cs  contributing  and  gives 

c  _j_  f  t-4-  («UI,t|/ci  -  s*“l,'yl/Cp)l  1  .  <«> 

■  -  -  4»p„‘  l  aiia’j  l  1  J 

In  some  applications  this  form  of  G  =  G,  +  G,  (i.e.  in  (x,  <o)  domain) 
is  desirable.  If  not,  the  inversion  (w  ->t)  proceeds  as  follows.  Referring 
to  (14)  we  have 


Gx(x.t;y) 


_  J_ 

2n 


.  -iwt  _  ,  . 

dto  e  G  (x,h>,y) 


00 

■sj  *• e"  [  ^7  ) 1 

—CD 

"  l!  Z |/=„>  I  • 


Similarly  inverting  GjU.wjy)  in  (16)  leads  to 


ok '‘f-s  S'S*:  •:  • 


du  e"iut[  -1- ■  — ill—  T-l_r  (e^U-yl/c,  -  ,i“U-l|/cp) 

L  4npw  3x  jdx j  - I 


_  1  B*  1  f  A  e  ib)t  T  i(i»|x-y|/c  iwjx-y|/c  ] 

-477707-  du— r  '-Zl  '  pJ 
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Recognizing  the  w  integral  as  the  second  integral  (in  t)  of  the  delta 


functions  6( t- |x-y  |/cs )  and  6( t-  |x-y  |/c_) .  we  have 


.  ,  „  %  1  a*  1  f 

=  4^37777  [-fit  1 


li-zl  li"Zl 

L(t  _  - >  -  L(t  -  - 

c  c 

P  8 


Zl  i 

r’JJ  • 


where,  if  H(t'-r)  is  the  Heaviside  function. 


L(t--c)  =  j  H(t'-r)  dt'  = 

—  CD 


0  ,  tit 


t-T,  t  >  r 


Noto  that  if  the  two  spatial  derivatives  were  to  be  carried  out  on  the  L 


functions  in  G,  one  would  obtain  the  delta  function  character  similar  to 
that  in  G1(x.t;y)  above.  Finally  G(x,t;y)  =  G4  +  G4  produces  the  free  space 


Green's  function  in  (x.t)  space. 
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THE  GREEN'S  FUNCTION  FOB  TUB  LAMB  PROBLEM 


The  traction-free  bonndary  conditions  at  the  x,  =  0  plane  can  be 
expressed 


Bn (xJ=0) 


M3 


0 

M3, 

P3, 

(X+2p)9 

-  0  , 


(17) 


x  =  0 
x 


where  dj  =  d/dxj.  We  are  now  seeking  solutions  to  equation  (3)  with  f  =  I 
for  Xj  1  0  and  satisfying  (17).  Because  of  the  restriction  (17),  we  cannot 
now  fully  transform  (x  — >  k)  the  problem.  Hence  we  begin  by  inserting  (k,  -> 
x,)  the  matrix  G(k,u;y)  given  by  (11)  and  (9);  this  will  provide  the  free- 
space  portion.  F.  of  the  Green's  function  we  seek  in  (k'.Xj.w)  space.  Then 
we  will  work  on  the  portion,  R,  defined  below,  that  will  cause  (17)  to  be 
satisfied;  and  our  final  Green's  function  will  be 

H  =  F  -  R. 

We  define  the  free-space  portion  in  (k' , x, )-space  by 


F( 


(D 

"  27  i  _  dk. 


(18) 


uj’«. 


1  T  . 

PP  + 


(c  k  -to  ) 

P 


— L__  UT  +  vvTll  . 

(c  |  -WV“  ■"  " 


Since  the  residue  theorem  is  to  be  applied  it  is  necessary  to  clearly 
identify  the  poles,  relative  to  k, ,  for  fixed  k’  =  (klfk,).  Consider  the 
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SLtJjUXs  ’/ulW 


i  «  v  tyt', 


imMmasMmmmmm 


(reciprocal  of)  first  tern  involving  a  pole  above. 


cV-V  =  c* (k*-  “_) 
P  P  c* 

P 


:*  (k* 


+  k' 


1 


c 

P 


=  c* (k*-  k*  ) 
P  *  »p 


where 


k 


*P 


sgn(w) 


(19) 


and  similarly. 


*.*  *  *,.*  .*  . 
ck-o>  =  c(k  -  k  ) 
s  s  j  >s 


k 


*s 


sgn  (<d) 


(20) 


The  sgn(o>)  will  be  explained  below.  In  this  notation  we  have 


F(k',xjw;y)  = 


(21) 


-ik' *y ' 
e  -  - 

2np 


CD 

1  “■ 


i  ,  a  i  , 
c  (k,  -  k,  ) 
P  *  »P 


T 

PP  + 


— i - |hhT  *  „T1 

,  L—  --  J 


c  (k  -  k  ) 

$  »  IS 


It  remains  to  select  the  proper  integration  paths  for  the  application  of  the 
residue  theorem.  There  are  four  poles,  k,  =  -  k,p  ►  -  kls  .  (Since  p  =  k/k 
and  etc.  for  h  and  v,  there  is  an  apparent  pole  for  k  =  k  +k,  =0.  A 
little  algebra  shows  that  there  is  no  pole  there).  Since  we  are  seeking  a 
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<»:v 


w 


8 


K.V  '  V 


l.,v» 


free-space  solution  due  to  a  point  source  (at  x  =  y)  we  want  our  solution 


F(k\x,.w;y)  to  be  'outgoing'  or  i,  -)  1  ••  (By  'outgoing'  we  mean  the 


following.  Should  our  F  be  inverse  transformed  (to  ->  t)  one  would  have  a 


composition  of  wave  forms  like  exp[-iwt  -  k,  (x,-y,))  ;  and  similarly  for 


k,s.  For  all  such  forms  to  be  outgoing,  for  t  >  0,  it  is  necessary  that 


sgn(*k,  (x,-y, ) )  =  sgn(w)  under  all  conditions).  With  the  specification 


sgn(k  )  =  sgn(k,  )  =  sgn(w)  in  (18)  and  (19),  it  remains  to  include  in  oar 
P  S 


integration  path  the  poles  in  (21)  by  the  rule: 


for  x  >  y  ,  include  poles  k  =  k  ,  k  j 

i  '  i  >  ip  is 


for  x  <  y  ,  include  poles  k  =  -k  ,  -k 

3  1  *  '  p  *  | 


With  the  integration  paths  determined,  the  residue  theorem  is  applied 


to  (21)  producing 


F(  k  ' ,  x ^  ,(i),  y)  = 


.  e_ *1  elk.plX»"y.l  ~~T  .  eik.slX._y.lr  rrl  .  -I  1 

1  2p —  — H - EE  + - ^ —  l-h-h  +  n  J 

cpk.p  csk.s 


where  p,  h,  ?  are  the  p,  h ,  y  vectors  with  k,  evaluated  at  the  proper  pole. 


That  is. 


P  =  ^  (k1'kt*k,p‘«n(x,-y,))T 

P 


h  =  h  =  i  (k1.-k1.0)T 


-tp-  <k1k,ss8n(x,_y,)  *  k*k,8s«n(xryi)  *  _lk'l  y 


where  ks  =  |(k1,k,,kl#)  |  and  kp  =  |(k1,k8,klp)  |  Note  that  since 


re. 


kJp  #  k,s,  the  2  *n<*  ▼  vectors  are  not  orthogonal.  However,  h  remains 
orthogonal  to  5  and  S'.  We  will  also  nse  the  notation 


p+  =  J_  (k  ,k  .  ±  k  )T 

-  kp  1  1  *P 

-1  =  kV  (±kik.,  •  ±kxk,s  '  -  lk'|,)T  • 


The  wave  vectors  with  the  +  are  associated  with  the  upward  moving  waves 

(*i  y  y*)  and  those  with  the  -  with  downward  moving  waves  (x,  <  y, ) . 

Next  we  apply  the  boundary  operator  B  to  F;  however,  we  do  this  in 

(k',x,)  -  space.  Hence  define  the  transform  (x'  — >  k’)  of  B  in  (17)  by 


B(k ' , x ^ )  = 


i.e.  the  operator  in  (17)  with  d,  ->  ik*  ,  d,  ik,- 

Before  applying  B  to  F,  we  will  decide  the  general  location  of  the 
source,  i.e.  y.  We  will  assume  that  y,  >  0  (source  above  the  x,  =  0  plane). 
In  the  common  y,  =  0  situation,  our  point  of  view  will  be  that  y,  >  0  first 
and  then  one  lets  y,  ->0.  Hence  as  we  apply  boundary  operator  B  to  F.  at  xt 
=  0,  we  use  the  expression  of  F  in  (23)  with  x,  <  y, .  This  results  in  p  = 
p~,  etc.  and  produces 


M3, 

0 

fMk, 

0 

m3, 

^k. 

ikkx 

ikk. 

(X+2y )d 

mVi 


mV**'',’  ■  ,  . 
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BP(*$-0)  =  B(k',xj)  F(k'  .x^  ,w;y) 


.  -ik'*v'  ~ik  (x  -y  ) 

B(k* ,x  )  »*  /  1  i - if  -1  ■  ‘  ,-p-T 

’  2»  [  cV  -- 

l  p  »p 


-ik  (x  -y  ) 
e  » s  »  i 

».» 
c  k 
s  *  s 


(hh1  +  v-v_T) 


2p  *.*  ip  K  * 

c  k 
p  *P 


+  -  f-1  B(k',-k,  )  (hhT  +  y-v-T) 

Csk.s 


where  B(k.-kJp)  is  B(k',x,)  with  the  d,  replaced  by  the  appropriate  value  of 
ik} ;  i.e.  -ik,p  in  the  first  case  in  (25);  similarly  for  B(k',-k,#). 

Since  F  accommodates  the  forcing  term  (8(x-y)5(t)I  in  (x.t)-space)  and 
has  BF ( x , =0 )  given  by  (25),  it  remains  to  find  matrix  R  satisfying  the 
homogeneous  equation  with  boundary  conditions  BR(x,=0)  =  BF(x,=0). 
Moreover,  as  discussed  above  for  F,  we  want  R  to  be  "outgoing"  as  x,  ->  +  ». 

The  homogeneous  vector  equation  in  (k',x,)  space  appears 


(L(k’,xj)  -  pul)  u  =  0 


where  L(k',x,)  is  L(k)  with  ik,  replaced  by  3,-  Equation  (26)  represents  a 
sixth  order  ordinary  differential  equation  in  that  it  consists  of  three 
second  order  equations  in  x,.  Its  general  solution  consists  of  linear 
combinations  of  six  vector  solutions;  however,  all  outgoing  solutions  (as 


'J  Vfr  *4<Y* v* 


*,  ->  ®)  consists  of  combinations  of  the  following  three  vector  solutions 


ik  x 

e  *P  *  p+ 


ik  x 
is  *  . 
e  h 


(27) 


ik  x 

=  e  s  *  v+ 


That  these  g^  solve  (26)  follows  from  earlier  work  on  F;  however,  it  is  easy 
to  verify  this  directly.  For  example,  putting  u  =  f  in  (26)  we  have 


ik  x 

(L(k',x  )-p«  I)  e  P  *  p+ 


ik  x 

e  *P  (L(k',kjp)  g+-pw*g+) 


=  e 


ik  x  s  +, 

ip  >  ( (X+2p)k  p  -pw  p  ) 


ik  x 

ip  *  ,  *.*  *.  +  „ 

i  r  p(ck-u)p  =0 

p  K 


where  we  used  equations  (7)  and  (20) i  and  in  the  last  two  lines 
k2  =  k  *+k*p.  The  key  step  in  this  verification  is  the  first  one,  in  which 
L(k ' , x, )exp(ik,p)  =  exp(ikIp)L(k’ ,klp) .  This  is  a  routine  calculation.  In 
the  same  manner  ga  and  g,  also  solve  (26).  The  other  three  solutions  are 
the  "incoming"  (for  x}  >  0)  solutions  obtained  by  using  p~ ,  v_ ,  -kJp  and  - 
k|t  in  (26).  Since  these  are  of  no  interest  to  us,  we  can  express  each 
column  of  Rlk'.x,;;)  »»  t  linear  combination  of  gt ,  g,  and  g,  .  Thus,  in 


R  ( k  '  ,Xj  ,«o;y)  =  4(x^;k')  C 


{U-U'U 1  c 


where  C  =  C(k',w;y)  is  a  "constant"  matrix  (relative  to  x3).  Our  remaining 
task  is  to  find  the  nine  elements  of  C  so  that  BR(x3=0)  =  BF(x3=0).  Of 
course  one  wants  to  avoid  solving  a  system  of  nine  equations  for  the  Cjj,  so 
we  proceed  with  some  care.  We  return  briefy  to  BF(x3=0)  i°  (25)  and  rewrite 


BF ( x  =0)  =  a  B'(-k  )p  P~^  +  a  B’(-k  )  f  hhT  +  v_v  T  | 

1  P  »p  -  -  s  i  s  L  —  --J 


where 


B '  (k  )  =  B(k  ' ,  k  )  . 

j  * 


-ik'*y’  ik  y 


2  2 
c  k 
p  >P 


-ik'-y'  ik,  y 
ie  e 

2  2 
c  k 
s  is 


Referring  to  (27)  and  (28)  we  have 


BR (x  =0)  =  ]  B ' (k  )p+,  B ' (k  >h,  B'(k  )v+  |  C 
j  [_  3  p  -  J  s  -  3  s  -  J 


We  now  anticipate  equating  BF(x3=0),  in  (29),  to  BR(x3=0)  work  on  matrix 


C.  Write 


or  p  +  a  (r.  h+r  v  , 
p  pp  -  s  hp  -  vp  -  ) 


_  t 

a  r  .  p  1  +  a  (r  h  +  r  vV 
"  T'h  £  s  hh  -  vh  -  ) 


or  p  ^  +  a  (r  h+r  v-)^ 
p  pv  -  s  hv  -  vv  - 


This  is  possible  since  p  ,  h,  v  are  independent  vectors.  (The  a  ,  as  were 


inserted  as  a  convenience,  i.e.  to  cancel  shortly  with  those  in  (29)  above.) 


Note  that  with  this  form  of  C  we  obtain,  using  (27)  and  (28), 


rik  x  ik  x  ik  x  1 

e  Jp3p+,e  *s  h.e  3s  v+  j  C 


ik  x 
.  3P  3 


p+for  p  +  a(r  h  +  r  v  )  l1 
-  L  P  PP  ~  s  hp  -  vp  -  J 


^  r  ■!>« 

+  e  ,s3h|arp-+a  (r  h  +  r  v”) 

-  L  P  ph  c  s  hh  -  vh  -  J 


s  hh  -  vh  - 


ilc  X 

+  e  *s  3  v+f  a  r  g"  +  a  (r  h+r  v“)  ]T  , 
-  L  p  pv  c  s  hv  -  vv  -  J 


where  ap  and  ag  are  defined  in  (30)  above.  This  follows  from  the  algebra 


fact  4.  above. 


It  remains  to  calculate  the  nine  reflection  coefficents.  The 


computational  advantage  to  the  form  in  (33)  is  that  it  expedites  the 


evaluation  of  the  unknown  r's.  We  now  refer  to  BR(x9=0)  in  (31)  where  C  is 


defined  in  (32)  and  equate  the  result  to  BF(x}-0)  in  (29).  In  particular. 


equating  coefficients  of  g  h^,  v  ^  leads  to  these  three  3  by  3  linear 


systems : 
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/.  .  r  .  4  «  4  m  •  « 

-  -  ,  -  •  <  A 


B'(-k  )p~  =  r  B'(k  )p+  +  B'(k  )  r  u  h  +  r  v+ 

*p  -  pp  tp  -  »s  t  ph  -  pv  -  J 


B'(-k  )h  =  r  B'(k  )p+  +  B '  (k  )  f  r..  h  +  r  v+  ] 

*  s  -  hp  ip-  » s  L  hh  -  bv-J 

B'(-k  )v~  =  r  B ' (k  )p+  +  B'(k  )  f  r  .  h  +  r  v+  ] 

JS  -  vp  3p  c  is  L  vh  -  vv  -  J 


(34a) 


(34b) 


(34c) 


Note  in  particular  that  the  three  systems  (34)  can  each  be  solved 
independently  for  the  unknowns  (e.g.  for  fpp*  rph*  rpv  *n  (34a)).  Also  note 
that  the  three  systems  have  the  same  3  by  3  matrix  of  coefficients. 


A  -  f  B’(k  )p+  .  B' (k  )h  ,  B' (k  )v+  1 
L  JP£  *  s  -  i  s  -  J 


This  fact  can  of  course  be  used  in  solving  for  the  r's  by,  e.g.,  inverting  A 
and  multiplying  the  left  sides  of  (34)  by  A-1.  The  terms  in  A  come  from 
evaluating  the  matrices  B'(-)  =  B(k',-)  and  vectors  p+  and  v+  at  the  values 
of  k  and  k  *  given  by  (19)  and  (20).  Similarly  the  quantities  in  the 

3  p  J  ® 

left  sides  of  (34)  are  evaluated  at  -k.  and  -k,  . 

*p  *  s 

When  the  three  systems  in  (34)  are  solved  we  obtain 


k  k'k 


r  =  -D'/D  , 
PP 


r  =  0  , 
hp 


4k  k  k 

*S  P  /v*  i 


ph 

0  , 

r  =  -4 
pv  k 

hh  = 

-1  , 

r  =  0  , 
hv 

vh 

0  . 

r  =  -D'/D  , 
vv 

(k|g-k'  )/D, 


where 


%  2  2  2 

D  =  4k  k  k’  +  (k  -k'  ) 
ip  •  s  is 


D'  =  4k  k  k'*  -  (k*  -k'V 
ip  ts  *s 
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Notice  in  (36)  the  various  aspects  of  deconpling  that  takes  place  in  the 
reflection  aatrix  R.  For  example,  the  middle  three  equations  in  (36)  state 
the  well  known  fact  that  an  incidence  h  li.e.  SB)  wave  gives  rise  (in  this 
setting)  only  to  a  reflected  h  wave.  And  similarly,  incident  p  and  v  waves 
only  produce  p  and  v  components. 

In  summary,  in  (k ' , x,  ,<fl)-space ,  the  Green's  function  for  the  Lamb 
problem  can  be  expressed 


H(k'  ,xt  ,e>;y)  =  F(k' ,xj  ,o);y)  -  R (k' , x*  ,w;y) 


(37) 


where  F  is  given  by  (23)  and  R  by  (33).  In  some  applications  this  form  of  H 
is  adequate,  while  on  other  occasions  it  is  preferred  to  have  H  in 
(x.t)space.  Unfortunately,  a  complete  closed  form  inversion  back  to  (x.t) 
has  not  been  accomplished.  Perhaps  the  most  successful  attempt  has  been  the 
ingenious  methods  due  to  Cagniard  and  deHoop  [5]  ,  in  which  the  inversion 
(k'  -*x')  can  be  reduced  to  a  single  (finite)  integral  representation  for 
H(x,t;y).  Illustrations  of  this  technique  can  be  found,  e.g.,  in  Aki  and 
Richards  [1] . 
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ABSTRACT 


This  paper  discusses  the  classical  Lamb  problem  for  the  elastic  wave 
equation.  The  motivation,  for  the  authors,  is  to  be  able  to  conveniently 
construct  Creen's  functions  (matrices)  for  later  use  in  formulating  and 
solving  various  inverse  problems.  For  example,  we  will  want  to  be  able  to 
solve  for  perturbations  from  constant  reference  densities  and  Lame 
parameters.  Hence,  Creen's  function  for  the  homogeneous  isotropic  equation 
is  discussed.  Due  to  the  scattering  taking  place  in  inverse  problems  it 
is  usually  impossible  to  retain  the  P-SV  and  SH  decoupling;  hence,  we  do 
not  pursue  this  decoupling  in  the  formation  of  the  Green's  functions  herein. 

While  nothing  conceptually  new  is  presented  here,  the  approach  is  a 
bit  different  and,  we  believe,  is  helpful  in  isolating  some  important 
issues.  The  approach  is  algebraic  in  nature  and  makes  heavy  use  of  several 
simple  facts  from  linear  algebra;  for  example,  the  spectral  decomposition 
of  special  matrices.  This  approach  facilitates  some  helpful  decoupling, 
particularly  in  solving  for  reflection  coefficients. 
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